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ABSTRACT 

Navier-Stokes calculations were carried out in 
order to predict the heat transfer rates on tur- 
bine blades. The calculations were performed us- 
ing TRAF2D which is a two-dimensional, explicit, fi- 
nite volume mass-averaged Navier-Stokes solver. Tur- 
bulence was modeled using Coakley’s q— u> and Chien s 
k — e two-equation models and the Baldwin- Lomax al- 
gebraic model. The model equations along with the 
flow equations were solved explicitly on a non-periodic 
C grid. Implicit residual smoothing (IRS) or a com- 
bination of multigrid technique and IRS was applied 
to enhance convergence rates. Calculations were per- 
formed to predict the Stanton number distributions on 
the first stage vane and blade row as well as the second 
stage vane row of the Rocketdyne Space Shuttle Main 
Engine (SSME) high pressure fuel turbine. The com- 
parison with the experimental results, although gener- 
ally favorable, serves to highlight the weaknesses of the 
turbulence models and the possible areas of improving 
these models for use in turbomachinery heat transfer 
calculations. 


NOMENCLATURE 

C Axial chord 

D Leading edge circle diam eter 

Fr Frossling number up / >/ Rep 

k Turbulent kinetic energy 

C Turbulent length scale 

M Mach number 

Nu Nusselt number 

p Pressure 
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Pr 

Prandtl number 

<7 

Variable in the q - u> equations ~\/k 

Re 

Reynolds number 

S 

Surface distance 

St 

Stanton number 

t 

Time 

T 

Temperature 

Tu 

Turbulence intensity 

u 

velocity vector 

V* 

Shear velocity 

y 

Distance from a solid surface 

y + 

Distance in wall coordinates, yv* /v 

Greek Symbols 

P 

Inlet angle 

7 

Specific heat ratio 

€ 

Turbulence dissipation rate 

6ij 

Kronecker delta 

K 

Thermal conductivity of the fluid 

P 

Viscosity 

V 

Kinematic viscosity 

U) 

e/lc 

p 

Density 


Subscripts 

0 Total condition 

1 Inlet condition 

D Diameter as characteristic length 

n Derivative normal to the surface 

2 Exit condition 

T Turbulent quantity 

id Surface conditions (w T a11) 

INTRODUCTION 

Accurate prediction of the airfoil heat transfer 
rates is critical in the design of modern gas turbines. 
Performing this task is difficult by virtue of the very 
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complex flow phenomena present in such turbines. Ef- 
fects such as turbulence or organized unsteadiness are 
modeled because they are impractical to directly sim- 
ulate due to limitations on computer memory and 
speed. Algebraic models usually work well for two- 
dimensional shear flows in local equilibrium and their 
implementation is easy and economical. However, in 
general, it may not be enough to prescribe turbulent 
scales by algebraic formulas. The time and length 
scales are dependent on transport effects and there- 
fore need to be found from transport equations. Two 
equation models strike a good balance between the 
complex multi-equation models which represent more 
of the physics of turbulence and the simpler models. 
The model equations in low Reynolds number mod- 
els are solved through the buffer layer to the wall. 
Low Reynolds number two-equation models are able 
to mimic transition induced by the free-stream turbu- 
lence. As the cases considered herein included tran- 
sition, this ability is investigated in this paper. The 
phenomenon of reversed transition or relaminarization 
can also be accounted for using these types of mod- 
els[l]. The features relating to transition and relami- 
narization are particularly useful in three-dimensional 
calculations where it is difficult to take advantage of 
empirical relations for the location of the start and 
extent of transition. The empirical relations are also 
often determined for two-dimensional flows which may 
limit their range of applicability. 

In the present work Coakley’s q - w [2] and 
Chien s [3] k — e low Reynolds number two-equation 
models are used. Previous applications [4,5,6] of the 
models to the prediction of blade heat transfer have 
been reported. The above models were chosen for their 
good numerical convergence properties. The model 
equations were incorporated in a two-dimensional cas- 
cade analysis code. For the rotating blade, the flow 
was solved in the relative reference frame and the 
body forces were neglected. The code utilized was 
the TRA F2D(TR A nsonic Flow 2D)[7] code which is 
a compressible, mass-averaged, Navier-Stokes solver 
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Shuttle Main Engine (SSME) high-pressure fuel tur- 
bine. Stanton number distributions were predicted 
for the first-stage vane and blade rows as well as the 
second-stage vane row. The results of the calculations 
were compared with the experiment al measurements of 
Dunn and Kim [9] which were obtained using a short 
duration measurement technique. 


DESCRIPTION OF THE ANALYSIS 

The flow was modeled using mass-weighted com- 
pressible Navier-Stokes equations and perfect gas 
equation of state. The variation of viscosity with tem- 
perature was assumed to follow Sutherland’s law[10]. 
The effect of turbulence was modeled using two low 
Reynolds number two-equation models. Solutions us- 
ing the Baldwin-Lomax algebraic model were also gen- 
era ted for comparison. 


Turbulence Models 


The turbulence was modeled using two-equation, 
low- Reynolds number turbulence models. Transition 
to turbulence is automatically mimicked by such mod- 
els. Chien s k — c and Coakley’s q — ui models were 
chosen for this work. The general formulation of the 
two models is written below following [2]: 

(psi)t + (psnij + q,j ) j = Hi 
V'J = ~(l‘ +f*t/Pri)sij , / = |,2 (1) 

For Chien ’s model the variables sj and so are: 

«i = k = q~ , s 2 = f = wk 


f'T = C /( Upk'/v 


The Reynolds stresses can be calculated from: 


' B . . 


D, E and F are damping functions 
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D= 1. - exp(— 0.01 15*/ + ), 


£ = 




F — exp{— 0.5,v + ) 


and 



also, 

Pr* = 1.0, Pr t = 1.3, C /4 = 0.09, C t = 1.35, Co = 1.8 


For the q — u; model, 

si = qr , so = u; 

/ir = CpDpk/w 

The source terms in the model equations are: 

//, = 1 [C„ DS/w 3 - |p/« - 1]/W« (4 ) 

//« = [Ci(c^/w 2 - |i»/«) - c 2 W 2 

D is a damping function defined as: 

D = 1 . — e.r;>(— o7?) 




Figure 1: Convergence history- q — u) model, 2 levels 
of multigrid 


Where Tu is the intensity of turbulence at the inlet 
of the computational domain and t/i is the inlet ve- 
locity. The value of e at the inlet is specified using 
the following relationship: 

c = u; k — k 3 t~/ ( 


where t is the turbulence length scale representing the 
size of the energy containing eddies. This length scale 
is usually not reported as a part of the experimental 
conditions. For cascade conditions the length scale is 
assumed to be equal to a certain percent of the pitch. 
For example Ilah [11] assumed a length scale equiva- 



No of Iterations 



Figure 2: Convergence history- q - u model, without 
multigrid 

terms were needed for the model equations. Variable 
coefficient implicit residual smoothing was adapted to 
the model equations by sett ing up the coefficients using 
the eigenvalues of the model equations. This allowed 
larger CFL numbers and therefore faster convergence 
rates. A CFL number of 5.0 was used for all the cal- 
culations except for k — c model's equations for which 
the CFL was in the range of 1.25-2.5. For the q — u 
model, it was possible to extend the multigrid capa- 
bility of the code to the model equations. Addition 
of second order dissipat ion terms was found to be un- 
necessary on the coarse grid for the model equations. 
Figures 1 and 2 show the convergence history for the 
q — u) model with and without multigrid. Figure 3 
shows the convergence history for a calculation per- 
formed using Chien’s model. The ordinate in those 
figures is the logarithm of the square root of the sum 
of the squares of the residuals over all the grid points 
divided by the number of grid points. The calculations 
were done for the first stage vane to be discussed later. 
The convergence history varies with varying values of 
exit Mach number and number of grid points among 
other factors. 


Computational Grid 

The discretization was performed on a non- 
periodic, C grid [7] generated by using a modified ver- 
sion of the Crape [14] code. With this type of grid t he 
number of grid points on t he suction surface of the air- 
foil is different than that on the pressure surface which 
leads to a much less skewed grid construction. 


CPU Requirements 

The CPU requirement per iteration per grid point 
was approximately 5.0E-5 second for the Baldwin- 
Lomax calculations with two levels of multigrid, 3.2E-5 


Figure 3: Convergence history- Jb — e model 

second for the q — w runs with one level of multigrid 
and 5.0E-5 second for the Chien’s model runs. The 
number of iterations needed was approximately 500, 
2500 and 4000 for the above models respectively. 


RESULTS AND DISCUSSION 

The Stanton number is defined as: 

Of _ *-^1 wall 

Pre/C p V re j(T w -To ) W 

where T w and Tq are the wall temperature and the inlet 
total temperature respectively. p re j, C P and V re j 
are the reference density, heat capacity and reference 
velocity. In the experiments the inlet to the first stage 
vane row is selected as reference. 

Given the two-stage turbine overall pressure ra- 
tio, the pressure ratio across each row of airfoils is 
needed as input to the N-S solver. The pressure ratio 
across the row as well as inlet total pressure and tem- 
perature and the airfoil Reynolds number were taken 
from the MTSB [15] flow analysis code as provided 
by Boyle [16]. The experimental runs were made for 
two values of total pressure, resulting in two values of 
Reynolds number for each airfoil. Calculations were 
made for both values of Reynolds numbers. Table I 
contains the flow conditions of the turbine airfoils con- 
sidered in this work. The free-stream turbulence levels 
for the three cascade rows are necessary for the calcula- 
tions. The turbulence level upstream of the first stage 
vane was experimentally determined but was guessed 
for the other two cascade rows. 

As stated earlier, the calculations were performed 
using the q — u and Chien’s k ~ c two equation models 
as well as Baldwin-Lomax algebraic model. For the al- 
gebraic model the location of transition to turbulence 
on the two sides of the airfoil was specified to best 
fit the experimental results. It is worth mentioning 
that there are a host of empirical relat ions that can be 
used [17] to predict the start and the length of transi- 
tion for use with algebraic models in two dimensions, 
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Figure 4: First stage vane and a typical grid 
however, that was not done here. 

First Stage Vane 



Figure 5: Variation of the grid spacing from the wall 
in terms of y + 



Figure 6: Surface pressure distribution, first stage vane 



Figure 4 shows the geometry of the first stage vane Figure 7: Variation of St with length scale, q - 

and the C grid generated. The mesh has 388x64 grid w mo del, Re . , = 3.6E5 
cells. The grid points are distributed such that there 
are 64 grid cells on the wake line on the pressure side 
of the airfoil and 32 grid cells on the wake line on the *>.„ 
suction side. Figure 5 shows the distance of the first 
point away from the wall in terms of the dimensionless u# 
parameter y + on the first vane for the lower Reynolds 
number run. As can be seen, this value is consistently Stnooo 
below unity. For the higher Reynolds number run the 
wall grid spacing was halved to keep the same y + dis- 
tribution. Figure 6 presents the surface pressure dis- so 
tribution on the blade surface as a function of surface 
distance. On the suction side there is an initial region M 
of favorable pressure gradient. The flow is likely to be 
laminar in this region for small values of free-stream 



turbulence. On the pressure side the flow encounters 
a small region of adverse pressure gradient following 
a short region of favorable pressure gradient. The ad- 
verse pressure gradient again can trigger transition for 


Figure 8: Variation of St with length scale, k 

e model, Re 2 = 3.6F5 
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Figure 11: First stage rotor blade and a typical grid 



Figure 12: First stage rotor, surface pressure distribu- 
tion 



s/c 


First Stage Blade 

Figure 11 presents the rotor blade of the first stage 
of the SSME turbine at the midspan. The mesh has 
388x64 grid cells. The grid points are distributed such 
that there are 64 grid cells on the wake line on the pres- 
sure side of the airfoil and 32 grid cells on the wake 
line on suction side. The same care regarding the grid 
resolution close to the wall was exercised when con- 
structing the grid for this blade. Figure 12 is a plot of 
the surface pressure distribution as a function of the 
surface distance. Figures 13 and 14 compare the Stan- 
ton number results with the experimental results. The 
laminar flow calculation for the high Reynolds number 
separated on the suction side. The Frossling analysis 
yields values of 0.019 and 0.014 for the low and high 
Reynolds number stagnation point Stanton numbers, 
respectively. The corresponding experimental values 
are 0.016 and 0.015 resulting in Frossling numbers of 
0.85 and 1.07. The location of transition on the suc- 
tion surface is further upstream of the experimentally 
determined location. Referring to figures 12-14, both 
of the models predict transition on the suction side 
through the first adverse pressure gradient. The suc- 
tion side level of Stanton number in the fully turbu- 
lent regime is correctly predicted for the low and the 


Figure 13: First stage rotor, surface Stanton number, 
low Reynolds number 



0 . 


S/C 

Figure 14: First stage rotor, surface Stanton number, 
high Reynolds number 
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Figure 16: Second stage vane surface pressure 


Figure 15: Second stage vane and a typical grid 

high Reynolds number cases with both of the mod- 
els. On the pressure side, the q - u model underpre- 
dicts the experimental Stanton number distribution. 
Chien’s model and Baldwin-Lomax model both under- 
predict the pressure surface heat transfer for the high 
Reynolds number condition. Dunn et al. [19] numer- 
ically investigated the effect of surface roughness on 
the heat transfer characteristics of this blade since the 
blade surface was rough. They concluded that the sur- 
face roughness effect was limited to the high Reynolds 
number condition where there was a small effect on 

the pressure side and a 12% increase on the suction 
side. 


Second Stage Vane 

Figure 15 presents the geometry and a typical grid 


ICO 



Figure 17: Second stage vane surface Stanton number 
low Reynolds number f 


*0 0 




Reynolds number cases were calculated t.o be 0.01.55 
and 0.0096, respectively. The experimental measure- 
ments are smaller than the above results and yield 
Frossling numbers of 0.75 and 0.9 which are quite low 
in the highly turbulent environment of the vane. The 
transition process has been captured with the q ^ 
model, albeit not as well for the higher Reynolds num- 
ber case. The k - ( model damped out the transition 
process, as was the case in the first vane. The data 
shows a phenomenon resembling a second transition 
process at the back of the blade on the suction sur- 
face which was not captured. The pressure side heat 
transfer results are satisfactory using both of t he two- 
equation models although Chien's model s agreement 
with the experimental measurements is better. 


SUMMARY AND CONCLUSIONS 

Navier-Stokes calculations to obtain the Stanton 
number distribution on the airfoil surfaces of Rocket- 
dyne SSME fuel high pressure turbine were performed. 
Coakley’s q - u> and Chien's k - < low Reynolds num- 
ber turbulence models and the Baldwin-lomax alg< - 
braic models were used. 1 he model equations wen 
solved using an explicit scheme. It was possible to ap- 
ply the multigrid technique to the solution ol the </-~ 
model. Implicit, residual smoothing was adapted to the 
solution of the equations for both of the two-equation 
Q f •di'we strategies helped to acccler- 
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